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SPECTRAL  MEASUREMENTS  OF  THE  WALL  SHEAR  STRESS  AND  WALL 
PRESSURE  IN  A  TURBULENT  BOUNDARY  LAYER:  THEORY 

INTRODUCTION 

Turbulent  boundary  layer  wall  pressure  fluctuations  have  received 
considerable  attention  over  the  past  30  years.  Willmarth  and  Wooldridge 
[1,2]  presented  spectral  and  correlation  measurements  of  the  flucuating  wall 
pressure  made  in  a  thick  turbulent  boundary  layer  of  air.  The  physical 
interpretations  made  by  Willmarth  and  Wooldridge  provided  much  of  the 
framework  on  which  subsequent  models  of  the  convected  wall  pressue  field 
were  developed.  The  issues  of  acoustic  contamination  (associated  with  the 
flow  facility)  at  low  frequencies  and  attenuations  at  high  frequencies 
resulting  from  spatial  averaging  (associated  with  the  finite  size  of  the 
face  of  the  pressure  transducer),  discussed  by  Willmarth  and  Wooldridge, 
continue  to  pose  severe  constraints  on  similar  experimental  research 
efforts.  More  recent  efforts  by  Bull  and  Thomas  [3]  and  Schewe  [4]  have 
been  aimed  at  resolving  the  pressure  flucuations  associated  with  the  small 
scale  (high  wavenumber)  structures  in  the  near  wall  region. 

Comparatively  few  investigations  of  wall  pressure  fluctuations 
performed  in  hydrodynamic  facilities  have  been  reported.  Carey  et  al.  [5] 
made  wall  pressure  measurements  in  a  turbulent  pipe  flow  of  water.  The  power 
spectra  were  shown  to  agree  with  the  measurements  of  Bakewell  et  al.  [6], 
which  were  made  in  a  turbulent  pipe  flow  of  air.  In  comparing  the  spectra, 
an  outer  flow  variable  nondimensionalization  was  used.  These  measurements 
were  also  compared  using  the  method  proposed  by  Corcos  [7]  to  correct  for 
the  effects  of  spatial  averaging  and  were  shown  to  agree  very  well.  In 
principal,  results  from  investigations  performed  in  air  may  be  scaled  to 
predict  the  spectral  levels  existing  for  a  particular  hydrodynamic  appli¬ 
cation,  if  corrections  for  spatial  averaging  are  taken  into  account  and 
changes  in  the  turbulent  structure  with  Reynolds  numbers  are  negligible.  To 
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investigate  boundary  layer  control  schemes,  such  as  polymer  additives, 
compliant  coatings,  or  large  eddy  breakup  devices,  for  reduction  of  flow 
noise,  the  use  of  water  as  a  working  fluid  is  necessary.  The  class  of 
hydrodynamic  boundary  layers  of  interest  have  sufficiently  small  viscous 
length  scales  to  cause  spatial  averaging  to  always  pose  a  limitation  in  the 
highest  frequencies  (or  smallest  scales)  which  can  be  resolved  with 
contemporary  transducers.  The  correction  method  proposed  by  Corcos  [7]  is 
approximate  in  nature,  as  discussed  by  Willmarth  and  Roos  [8].  The  method 
does  provide  a  systematic  correction  procedure  that  is  useful  for  comparing 
results  from  various  investigations.  In  addition,  for  a  particular  boundary 
layer  and  transducer,  the  method  may  be  used  to  estimate  the  highest 
frequency  for  which  an  acceptable  signal  to  noise  ratio  exists. 

The  related  problem  involving  the  flucuating  wall  shear  stress  beneath 
turbulent  boundary  layers  has  received  considerably  less  attention.  A 
recent  discussion  of  the  current  knowledge  of  wall  shear  fluctuations  is 
given  by  Alfredsson  et  al .  [9].  Measurements  of  the  ratio 
have  varied  drastically  between  various  investigations.  The  spectra  of  the 
wall  shear  stress  flucuations  have  not  been  well  established.  The  inherent 
difficulties  associated  with  the  various  measurement  methods  are  largely 
responsible.  Alfredsson  et  al .  concluded  that  flush  mounted  hot  film  sensors 
may  be  used  to  measure  the  mean  and  fluctuating  wall  shear  stress  in  water 
applications  where  the  heat  loss  from  the  hot  film  to  the  substrate  is  not 
significant.  The  validity  of  this  measurement  method  requires  that  the 
thermal  boundary  layer  developed  over  the  hot  film  be  entirely  within  the 
viscous  sublayer  of  the  turbulent  boundary  layer.  In  addition,  the  unsteady 
thermal  response  of  the  sensor  poses  a  constraint  on  the  highest  frequency 
wall  shear  fluctuations  which  may  be  accurately  resolved.  This  report 
contains  two  derivations  necessary  in  order  to  determine  the  applicability 
of  a  particular  hot  film  sensor  for  a  particular  application.  The  first 
derivation  involves  a  steady  state  heat  transfer  analysis,  from  which  the 
thermal  boundary  layer  thickness  is  determined.  The  second  derivation 
involves  an  unsteady  heat  transfer  analysis,  from  which  the  one  half  power 
frequencies  are  determined. 
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The  experimental  results  from  this  investigation  are  given  by  Keith 
and  Bennett  [10],  The  goal  was  to  resolve  both  the  wall  pressure  and  wall 
shear  stress  spectra  in  a  fully  developed  hydrodynamic  turbulent  boundary 
layer.  In  addition  to  the  frequency  response  of  the  hot  film  sensors,  their 
absolute  calibaration  posed  an  additional  constraint  on  the  accuracy  of  the 
measured  spectra.  This  issue  is  discussed  in  detail  in  [10],  With  regard 
to  the  results  presented  here,  it  is  emphasized  that  for  hydrodynamic 
applications,  both  the  thermal  boundary  layer  thickness  and  the  frequency 
response  associated  with  flush  mounted  hot  films  must  be  considered  at  the 
outset  of  any  experiment.  For  many  applications  in  air,  the  thermal 
boundary  layer  thickness  is  an  order  of  magnitude  less  than  the  viscous 
sublayer  thickness.  For  the  hydrodynamic  application  here,  the  thermal 
boundary  layer  was  approximately  one  half  as  thick  as  the  viscour  sublayer. 
Attenuations  resulting  from  the  thermal  response  limit  the  frequencies  for 
which  the  spectral  levels  are  accurately  resolved. 

Currently,  it  is  not  possible  to  estimate  the  effects  of  spatial 
averaging,  due  to  the  lack  of  measurements  of  the  cross  correlation  of  the 
wall  shear  stress  field.  For  the  boundary  layer  investigated  here,  the 
streamwise  length  of  the  hot  film  sensors  was  an  order  of  magnitude  smaller 
than  the  diameter  of  the  pressure  transducers.  As  a  result,  the  size  of  the 
pressure  transducers  was  the  limiting  factor  with  regard  to  the  frequency 
range  over  which  both  the  wall  pressure  and  wall  shear  stress  were 
accurately  resolved.  Significantly  more  experimental  measurements  of  the 
wall  shear  stress  field  are  needed  in  order  for  an  understanding  comparable 
to  the  wall  pressure  field  to  exist. 


3 


TR  8295 


DERIVATION  OF  THE  STEADY  STATE  RESPONSE  OF 
THE  HOT  FILM  SENSORS 


Nomenclature 


MX) 

r“  e  U^e/v 
ii(y) 

T 

OD 

T(y) 

T^(y)  =  T(y)  -  T^ 
k 

S 

p 

K  =  k/pCp 


Turbulent  boundary  layer  thickness  (in.)t 
Free  stream  velocity  (ft/sec), 

Momentum  thickness  Reynolds  number. 

Mean  streamwise  velocity  (ft/sec). 

Thermal  boundary  layer  thickness  (in.). 

where  ^t^’^o  =  *tL’ 

Fluid  temperature  (*F)  for  y  >  4^(x), 

Mean  fluid  temperature  (®F), 

Thermal  conductivity  (BTU/(sec-ft^  -  “F/ft)), 
Constant  pressure  specific  heat  (BTU/(slug  -  “F)), 
Fluid  density  (slug/ft^). 

Fluid  thermal  diffusivity  (ft^/sec), 

2 

Heat  transfer  rate  per  unit  area  (BTU/(sec  -  ft  )), 
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2 

V  Dynamic  viscosity  (Ibf  -  sec/ft  ), 

2 

V  Kinematic  viscosity  (ft  /sec), 

2 

Tq  Mean  wall  shear  stress  (Ibf /ft  ), 

0  =  Cpv/k  Prandtl  Number 

u*  =  (Tq//»)^^^  Friction  velocity  (ft/sec), 

L  Streamwise  length  of  the  hot  film  (in.)- 

Following  Bellhouse  and  Schultz  [11],  the  thermal  energy  integral 
equation  may  be  expressed  as 

u(y)T„(y)<ly  =  ,„(x)/pCp  =  -  K  ^ 

From  [11],  a  steady  thermal  boundary  layer  profile  of  the  form 
To(n)  =  b(l-n)^  (l+n)  =  b(l-n)^  dV) 
is  used,  where 

b  =  T(o)  -  T  =  constant, 
and 

n  »=  y/«^(x). 

Then 


^^0  an  _  b 

an  ay  “  a^(x) 


[-2(l-n)(l-n^) 


(l-n)^(-2n)] 


and 


-2b 


(2) 
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Next,  it  is  assumed  that  the  thermal  boundary  layer  thickness  6^(x)  is 
entirely  within  the  viscous  sublayer  of  the  turbulent  boundary  layer,  such 


u(y)  =  y/y  . 


Combining  Eqs.  (2)  and  (3)  with  Eq.  (1)  results  in 


a  ^t^^^  2  2 

f-  r  T  vb(l-n)‘'(lV)dy 

ax  0  0 


fr  I  —  C4?{x)(l-n)^(l-n^n3dn. 
ax  0  y  X 


Assuming  xq  to  be  constant  over  the  region  Xq  <  x  <  Xq  +  L  gives 


which  upon  integration  yields 


T^f-r  =  2«,(x)  d  ‘t'*’ 

V**  ‘ 


Integrating  Eq.  (5)  from  xq  to  xq  +  L  yields 


Xa'*'  L  X 

0  T  ^+*1  O 

I  K  dx  =  :f|-  4^(x)d4^(x), 

Xq  0 


which  results  in 


KL  »  Tq  ^5y 
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or 

=  45yK/(T^L2). 


(6) 


Defining  =  L  u*/»,  Eq.  (6)  may  be  expressed  in  the  form 


^  /I  /iic/  il/3,,+v-2/3 

=  (45/o)  (L  ) 


(7) 


In  this  particular  investigation  [10],  L  =  0.005  (in.)  and  a  =  6.40 
(corresponding  to  a  fluid  temperature  of  72*F).  The  thermal  boundary  layer 
thickness  was  entirely  within  the  viscous  sublayer  for  each  Reynolds  number 
investigated,  as  shown  in  table  1.  Here,  the  viscous  sublayer  thickness  is 
assumed  to  be  10  viscous  lengths. 

Table  1.  Boundary  Layer  and  Hot  Film  Parameters 


U«  (ft/sec) 

Re 

6  (in.) 

i>/u*  (in.) 

Lu*/ V 

AtL^^*/ w 

10 

8,200 

0.93 

0.00034 

15 

4.65 

20 

13,400 

0.74 

0.00019 

25 

5.72 

Having  established  that  Eq.  (3)  holds  throughout  the  thermal  boundary 
layer,  an  expression  relating  the  mean  wall  shear  stress  to  the  steady  state 
mean  heat  transfer  may  now  be  derived.  The  steady  state  mean  heat  transfer 
rate  (with  respect  to  the  streamwise  element  length  L)  per  unit  area  is 
defined  as 


dw  •  L  I  qvrf^x)dx.  (8) 

*0 

From  Eqs.  (1 )  and  (2) , 

=  2IC  b/a^(x)  , 
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and  from  Eq.  (5) , 


pCp  “  15w  dx 


Combining  Eqs.  (,9)  and  (8)  yields 


_  2pC  T  b  a 

ow'i^r  ‘tO'X'V’') 


*tl 


From  Eq-  (6) , 


pC  T  b 

—  DO 


45uKL 


Eq.  (10)  may  be  expressed  in  the  form 

^  =  k  b  (3/(5vLK))^^^  (11) 

which  is  the  desired  relationship  between  qv^  and  tq,  Eq.  (11)  relates 
the  steady  state  mean  heat  transfer  rate  from  the  hot  film  to  the  mean 
wall  shear  stress  t^.  In  practice,  it  is  not  the  mean  heat  transfer 
rate  that  is  determined  in  order  to  calibrate  the  sensor,  but  rather  the 
anemometer  output  voltage  E  necessary  to  control  the  hot  film  at  a  constant 
temperature.  Sandborn  [12]  has  shown  that  the  relationship  between  the  mean 
wall  shear  stress  and  anemometer  output  voltage  is  of  the  form 

=  AE^  +  8.  (12) 


where  E  is  the  voltage,  and  A  and  8  are  calibration  constants.  While  the 
accurate  determination  of  these  constants  for  particular  applications  re¬ 
quires  considerable  effort,  the  result  yields  an  equation  relating  the 
measurable  quantity  E  to  the  desired  quantity  t^.  The  particular  method 
of  calibration  used  in  this  investigation  is  fully  discussed  in  [10]. 
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DERIVATION  OF  THE  FREQUENCY  RESPONSE 
OF  THE  HOT  FILM  SENSORS 


To  determine  the  dynamic  response  of  the  hot  film  sensors,  the 
quantity  T^(y)  +  T^(y)e^“^  is  combined  with  the  unsteady  form  of  the 
thermal  energy  integral  equation  to  obtain 


a_ 

at 


r 

0 


(Tj,(y)  +  T^(y)e^"^)dy 


a  4„-{- 

a7  ^  Tf(y)e’“'')dy 

•  ,.o- 

The  first  term  in  Eq.  (13)  may  be  reduced  to 

iait  ^t^^^ 

i«  e^“^  T,(y)dy. 

0  ^ 

Following  [11],  the  fluctuating  thermal  boundary  layer  profile  is 
approximated  as 

T^(n)  =  -b[B(«)n(l-n)^l-*-2n)]. 

where  n  =  y/a^(x).  The  steady  thermal  profile  is  approximated  as 
TgCn)  =  b(l-n)^(l-n^). 

and  the  fluctuating  streamwise  velocity  u(n)  is  expressed  as 


(13) 


(14) 


(15) 


(16) 
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u(n)  =  (Tq/v  +  e^“^/y)Ti  4^(x).  (17) 

Combining  Eqs.  (15),  (16),  and  (17)  with  Eq.  (13)  and  performing  the 
necessary  integrations  results  in 

-  i«  B(«)4t(x)  (3/20)  e^“t 

+  CV^  ^  [fg  4^(x)  (1-B(*.)e^“^)] 

=  (2+8(0.)  e^“^).  (18) 

Equating  the  constant  terms  in  Eq.  (18)  to  zero  and  integrating  results  in 
4^(x)^  =  45KyX/TQ,  (19) 

which  is  consistent  with  Eq.  (6)  from  the  steady  state  analysis. 

4  ▲ 

Equating  the  coefficients  of  e  in  Eq.  (18)  to  zero  and  neglecting 
higher  order  terms  results  in 

n  n  44. (X)  -  _ 

15  *t(^^  aT"  ■  20 

Eq.  (20)  may  be  integrated  over  the  sensor  length  L  (with  respect  to  the 
streamwise  direction  x)  to  obtain 

B(«)  =  r  +  40  (21) 

0 

where  the  relationship 

,  34.(x) 

<<■»  « 
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obtained  In  the  steady  analysis  has  been  used,  and 


Combining  Eqs.  (21)  and  (15)  and  defining  qi(x)  as 
3T^ 

y«0 


results  In  an  expression  for  the  fluctuating  heat  transfer  given  by 


„  r3  ^  31“ 

~  6^(x)  ^2  40K  *t^ 


(22) 


The  mean  value  of  qi(x)  defined  as 


’l  =  L  > 


is  obtained  from  Eq.  (22)  as 


;r  rh,  3i«  .2  ,-i  dx 

‘’l  ■  L  Tq  40K  *t  ^  J  «^(x)- 


-:r  ,  L 


(23) 


The  mean  value  of  the  steady  state  heat  transfer  q^  may  be  expressed  as 


—  2^  f  dx 

^0  “  L  J  i^(x)  • 


(24) 


Combining  Eqs.  (24)  with  (23)  yields 


TT/TT  !l-  rl  *  ^  x2  i-l 

■  2t„  *^2  *  40K  *t  ^ 

0 


(25) 
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Taking  the  magnitude  of  Eq.  (24)  yields 

0 

From  Eq.  (19), 

«^(x)  =  (K45vx/t^)^^\ 

Then 

l2  If  ,^Ac  V/  ^2/3Hv  3.K4SwL.^^^ 

it  =  r  J  (K45vx/t  )  dx  =  f(  ,  ) 

I  L  0  o  ^  ^0 

From  Eq.  (26),  the  one-half  power  frequency  expressed  as 

“1/2  ^t  /  (20  K)  =  1. 


(26) 


(27) 


(28) 


Combining  Eqs.  (27)  and  (28)  and  rearranging  the  terms  results  In 

-7  ,  ,  5  -1/3 

“1/2  =  (LvVu*‘=)‘=]  .  (29) 

Eq.  (29)  defines  the  one-half  power  frequency  for  the  hot  film  sensors. 

Table  2  presents  both  the  dimensional  and  nondimenslonal  one-half  power 
frequencies  for  this  particular  investigation. 


Table  2.  Hot  Film  Sensor  One-Half  Power  Frequencies 


U  (ft/sec) 

QB 


fl/2  “1/2 


“l/2  w/(u*)‘ 


10 

20 


468  3.18 

1016  2.40 


0.24 

0.16 
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CONCLUSIONS 

The  analysis  pertinent  for  the  use  of  flush  mounted  hot  film  sensors 
to  measure  the  fluctuating  wall  shear  stress  has  been  presented.  For  the 
case  where  the  thermal  boundary  layer  thickness  is  of  the  order  of  or  less 
than  the  fluid  viscous  sublayer  thickness,  the  mean  output  voltage  from  the 
sensor  is  directly  related  to  the  mean  wall  shear  stress.  An  approximate 
two-dimensional  analysis  of  the  unsteady  heat  transfer  results  in  an  ex¬ 
pression  that  characterizes  the  frequency  response  of  the  sensor.  Use  of 
this  equation  for  typical  applications  in  water,  i.e.,  those  requiring  a 
relatively  low  sensor  temperature,  shows  that  the  attenuation  at  high 
frequencies  due  to  the  thermal  response  is  significant.  The  results  of  this 
analysis  show  that  the  flush  mounted  hot  film  sensors  used  for  the  investi¬ 
gation  of  Keith  and  Bennett  [10]  were  effective  in  resolving  the  lower 
frequencies  of  the  fluctuating  wall  shear  stress. 
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